A Generalization of Heaviside's Expansion Theorem 

By W. O. PENNELL 

The expansion theorem is one of the most frequently used methods of 
evaluating operational forms arising from the operational calculus de- 
veloped by Heaviside. The original theorem, however, is applicable, in 
general, only to expressions containing integral powers of the operator d/dt. 
This paper describes an extension to, or a generalization of the original 
expansion theorem whereby, in general, operational forms with either 
fractional or integral powers of the operator can be evaluated. A number 
of operational equivalents are given to be used with the theorem, one of 
which is the equivalent used by Heaviside. Examples of the application 
of the theorem to electric circuit problems are shown. 

THE well known expansion theorem given by Heaviside in Vol. II 
of his "Electromagnetic Theory" may be stated as follows: 
An operational equation of the form /; = Y(p)/Z(p), may under 
certain well known restrictions on the functions Y and Z, have as its 
solution 

* = S+?^'- ■-"-- « 

p is the differential operator d/dt, and pi, p% • • • are the roots of 
Z{p) = 0. Z'(p n ) is the result of substituting p n for p in d(Z{p))/dp. 
The theorem is true only when no root is zero and all roots are unequal. 
Y(p) and Z(p) must contain p to positive integral powers only. Various 
proofs of this theorem have been given and perhaps the simplest de- 
pends upon the expansion of Y(p)/Z(p) by partial fractions. 

The expansion theorem is valuable in the solution by ODerational 
methods, of problems in mathematical physics, and especially electric 
circuit theory problems. 

Generalization of the Expansion Theorem 

The generalization of this theorem may be stated as follows: Under 
certain circumstances it may be possible to write the operational 
equation 

" Z{p) aS D(q)' 

where q is a function of the operator p. With suitable restrictions on 
the functional forms of N and D the solution of the operational equa- 
tion is given by 

_ 7V(0) N{g n ) 

k ~ D{0) + \q n D'{q n ) W ' qn)t {2) 

where if/(t, q„) is the equivalent of the operational expression q/(q — q n ) 
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and gi, q 2 , • • • represent the roots of D(q) = 0. If q is the differential 
operator, that is if 

/. d 

then as is well known 

-g— = P = 0* 

Q — Qn P — Pn 

and (2) becomes the Heaviside expression (1). 

A proof of the generalized theorem equation (2), is as follows: 
By a theorem of partial fractions: 

N(q)_ N( gi ) N(q 2 ) N(q n ) 

D{q) (q - qi)D\ qi ) "*" (q - qJD'fo) * " ' "*" (q - q n )D'{q n ) {6) 

where q u q 2 , • • • q n are the roots of D(q) = 0. The above theorem is 
true when D(q) and N(q) are rational polynomials and N(q) is of a 
lower degree than D(q). Further limitations are that no root can be 
zero and all roots must be unequal. 

In writing the above identity in terms of operators it is tacitly as- 
sumed that the operators obey the three fundamental laws of algebra, 
the associative, commutative and distributive laws. 

Now 

_J_=_i + «_. (4) 

5 - 2» ffn q n (q - q n ) 

Substituting (4) in (3) 

N(q) _ N( qi ) N(q 2 ) N(q n ) 

r 7 — _ n r.//. V ~r " " • -r 



D{q) (- qi )D'( qi ) {-q 2 )D'{q 2 ) ' ' {-q n )D'(q n ) 

n ZnD'{q n ) \q - q n 
_N(0) N(q n ) 

-m + ^^jyM+ (t ' qn) ' (6) 



where 



^= *«.*.>• 



The expression fails where N(0)/D(0) is infinite. When the operator 

^ d 
q=P = dt 
then 

P-Pn 

and (6) becomes the Heaviside Expansion theorem. 
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Although the above proof of (6) is for cases where D(q) is a polynom- 
ial, if D(q) is a transcendental function which can be expanded by the 
process shown, the equation will still hold. It is shown in treatises on 
trigonometry that tan at, cot at, 1/sin at, 1/cos at, 1/sinh at, 1/cosh at, 
tanh at, and coth at, all can be expanded in an infinite series of partial 
fractions which are identical l with the expansions obtained by applying 
the process of equation (3). 

Equivalents to be Used in Generalized Theorem 

In applying this theorem the following operational equivalents are 
useful : 
Equivalent No. 1: 

Let 

A d 

Then 

a _ p 



q — a p — a 



= e»'. (7) 



This is the equivalent used in the expansion theorem by Heaviside. 
Equivalent No. 2: 



Then 
where 



- ph '" = e a2l \_\ + erf (a/" 2 )] (8) 



Sl/2 _ 



q - a p 

,1/2 



Equivalent No. 3: 
Let 



rat 111 
erf K' 2 ) = -== I er»d\. 

q = pu» = I \ s = a positive integer. 



Then 



q-a p 
where 



1 Except in some cases for the first term F(0)/Z(0). 
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1 /Vi s " 

1 i , \>i 



— + 1 1 J ° 



r(- — - + 



Equivalent No. 4: 
Let 

Then 

— ^— = -fi— = cosh a 1 ' 2 /. (10) 

q — a p- — a 

Equivalent No. 5: 

Let 

.-,-(«'■ 

Then 

^— = (1/3)€» I/8 ' + (2/3)e-» ,/3 ''2 C os (a l W?j ) • (U) 



g — a p A — a 

Equivalent No. 6: 
Let 

-*-«)'■ 

Then 

g _ £ 4 _ i 
q — a p 4 — a 



„ = | cosh (a 1 / 4 /) + § cos (a 1 ' 4 /). (12) 



Equivalent No. 7: 
Let 

d . , \ 1/2 



Then 



j -<* + *)"-(£+ ft)' 



2 - a 6 - a 2 6 - a 2 * + b - a 2 * 2 " 1 ^ w; 



- 1 -^e^-«'erf (aVF). 
o — a 1 



where b p^ a 2 

Equivalent No. 8: 
Let 



2 = 



£ + 6 



486 BELL SYSTEM TECHNICAL JOURNAL 

Then 

_J_ = —L- \ **«/a--> + ae-wio (~) 
q — a 1 — or L \ 2 / 

+ ^ 2) 2 e° 26t/(l - a2) f" e<«»+i>w/2(«*-i)j / *j \ ^ 1 
where a 2 5^ 1 



and Jo ( -r ) = Jo ( -y ) = Bessel Function of the first kind. 

The above equivalents can be obtained by known operational 
methods and their derivation will not be given here. 

In the application of the generalized theorem to electrical problems, 
equivalents No. 1, No. 2 and No. 3, especially No. 1 and No. 2, are the 
ones which will be most frequently used. Equivalents No. 4, No. 5, 
and No. 6, since they involve only integral powers of p are of use in 
reducing the labor of applying the original expansion theorem to ex- 
pressions containing only these powers of the operator p or multiples of 
these powers. Their use in such cases is illustrated by example No. 3 
below. 

Equivalent No. 7 enables expressions like the following to be evalu- 
ated in closed form. 

1 1 



p + c(p + b) 1 ' 2 + d' {p + b) 3 ' 2 + cp + d' 

(P + &) 1 ' 2 1 

(p + b) 3 ' 2 + cp + d' cosh {p + by> 2 ' 



etc. 



In applying equivalents No. 2 and No. 7 some of the following proper- 
ties of the error function are often conveniently used. 

erf (-/) = - erf (/) and erf (*/) = 4= f e * d ^ 

■VTT Jo 

also 

|erf[*(0] = ^-^> ] V'(0 

and 

^.erf (at 1 '*) = ^'ar 1 ' 2 . 

at -yjjj- 

The value of erf (/) for different values of / may be obtained from tables 
of the probability integral as for example Pierce Table of Integrals. 
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The values of erf (it) for values of t from .01 to 2 are given in a table in 
London Mathematical Society Vol. 29, 1897-98, page 519. The values 
of erf (te iTli ) and erf (te~ ir,A ) are given by the following formulae: 



where 

and 

and 
and 



Tables of the values of these two integrals known as the Fresnel 
Integrals are given in various handbooks such as Jahnke and Emde. 

Examples of Application of Theorem 

A few applications of the theorem will be given. 
Example 1: 

The operational solution for the current entering an infinitely long 
ideal cable with a given impressed voltage of the form Ee~ at is 

£3/2 

K P 



erf (Je iTl4 ) = j2i[C(H2/ir) - iS(M2/ir)l, (14) 

erf (te-*"*) = V2l[- iC{tjIfn) + S{tJ2fr)'], (15) 

Xt^fif* / _/2 \ 

sin (^)dt, 

C(- ihj2/r) = - *C(/V2/ir) 
Si-it-fifr) = iS(t^2~fir). 



P +a' 
K being a constant and 

A d 

To evaluate p v<l l(p + a) in closed form call p u - = q then 

£3/2 ^ q s 
p + a q- + a 

Since the theorem applies in general only when the degree of the 
numerator is less than that of the denominator we will write 

q 3 aq 

= q - 



q' + a q 2 + a 
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and 

g _ Y(g) _ F(0) Y(g n ) 

q 2 + a Z{q) Z(0) ^ f <Z«Z'(<Z„) n? "' A 

TO - n y (g») - g" - 1 

Z(0) U - q n Z'(q n ) 2q 2 2q n ' 

Z(q) = q 2 + a and the roots of Z(q) = are q n = ± ia JI2 , 

^q ni t) = e qnH [\ 4- erf ?„/ I/2 ] (see equivalent No. 2). 
So 

_^_ = _J— e -a<[i + erf (m 1 ' 2 / 1 ' 2 )] - T 4 T75 <r»'[l + erf (- m 1 ' 2 / 1 ' 2 )] 
q 2 -\- a 2ta 1 ' 1 lta vli 



= -Aj5<r at er{(ia l iH 1 ' 2 ). 



Hence 



q 2 + a s q 2 + a Y{Yi) i 

since 

AL/I / d Y" rl/2 

Example 2: 

The operational expression for the current entering at time 2 in a 
cable of distributed resistance R and capacity C with an electromotive 
force sin cot impressed is given by 

. cop 3 ' 2 

where 

*- d 

Put q = p 1 ' 2 . Then 

£3/2 q 3 



p 2 + CO 2 S 4 + CO 2 

Here 



TO 



-n ^n) _ gn 3 _ 1.. 



Z(0) "' ffnZ'foO 4g„ 4 4g„ 

If 

gA + w 2 = o, g „ = w l 'V ( - /4 >, tfi*e-*W* t - w 1 ¥ T,4) l - tfr** 4 *. 
So 
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^ t^- 2 = « VCT^ I^e""[l + erf (q H fti*)l 

= W ^ [ I^Mi g<a " I * + erf (« l/2 ^ /4> ^ /2 ) J 

" 4^W) '*" I ! + erf C- wlWfll » 



CO 



1/2 



_J_ g -«Mf-»/4) er f ( a ,l/2g-i(W4)/)/2)] 



= J^5[sin (c0t)S{co"H l l^2lir) 

+ cos («0C(tt"W B V2^r)]. 

The last transformation is obtained by means of formulae 14 and 15. 
Example 3: 
Evaluate 

= 1 

y £ 4 - 3£ 2 + 2 ' 

This can be solved by the expansion theorem in the usual way. A 
somewhat shorter method is to use the generalized theorem with the 
operator a = p 2 . Then 

g 2 — 3g + 2 * n 2g„ 2 - 3q n 
q n = 1,2; iK ffBI = cosh^ 1 ' 2 /. 

See equivalent No. 4. So 

y = \ + | cosh /V2 — cosh L 

Example 4: 
Evaluate 

sinh^ = sinh&ff = b sinh bq n erf /t 

sinha^ 1 ' 2 sinh ag a ag n cosh aq n 

— sinh aq= i sin tag. 
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The roots 1 of sin iaq = are 

flTT 

2„ = — , »=±1,±2, etc. 

id 

Substituting these values of q n in above we get 

. , mrb 

smhb 1= b + -5" r *~+*\ 1+ea (¥**)] 

n = :t2— cosh — 
etc.. 1 1 

A ? (-l) n sin— - 

— _ _|_ Z V _g-(n*ir»/o«)l. 

« Tn=1.2.3. ••• » 

If (sinh 6/> l/2 )/(sinh a/> l/2 ) is solved by the expansion theorem and the 
summation is extended over both positive and negative roots, the result 
is 

. mrb 

J A Sin 

- + - £ (- 1)« — e-f" 1 "/" 1 )'. 

7T „=!, 2. 3. • • • « 

In other words the summation quantity is just double what it should be. 
In order to correct this in practice, those who have used the theorem 
for such cases have extended the summation only over the positive 
roots, notwithstanding the fact that in similar cases with integral 
exponents such as, for example, 1/cosh ap the summation is extended 
over all the roots. The truth is the original expansion theorem is not 
applicable if either numerator or denominator contains p to a fractional 
form. In the above case were the problems to evaluate (sinh bp 213 / 
sinh ap 2 ' 3 ) the expansion theorem gives an entirely incorrect answer, 
while the correct answer is obtained from the extension to the theorem. 



Example 5. 



Here 






[i pen* 111 1 /»«„»»'» 

?»=±1, 

1 Excluding the root q n = which is not used in this case. 



HEAVISIDE'S EXPANSION THEOREM 491 

1 ! ^'Jo 



- * <r x vx. 



Example 6: 

If the problem is to evaluate 

sinh fr(j> + c) l/2 
sinh a (/> + c) 1/2 

Equivalent No. 7 is used. The details will not be worked out since 
they are quite similar to Example No. 4. The answer is 

sinh6(£ + r) 1/2 _ b 
sinh a(p + c) 112 a 

nirb 

2 2 



_ sin , 

+ 2 - Z (-!)"■ " 



7r„ = l,2,3... W 



W 2 7T 2 fl 2 C + n 2 ir 2 
C " -,2 



fl- 



If c = the above equivalent reduces to the answer of Example 

No. 4. 

Final Remarks 

Operational methods were used by Euler and other mathematicians 
prior to Heaviside. Their use, however, depended in general upon a 
formal definition of the operator. Heaviside, on the other hand, 
adopted a different procedure. In the differential equation of the 
problem he replaced the operator d/dt by p and obtained the solution of 
the resulting algebraic equation. He then determined the significance 
of the operator by the condition that it should give the complete solu- 
tion of the original differential equation subiect to equilibrium bound- 
ary condition. 

While Heaviside developed the operational calculus in a fairly 
workable and complete form he failed to correlate it or reconcile it 
with conventional mathematics or to put its theorems on a rigorous 
basis. The development since Heaviside's day has been due to a 
considerable extent to the engineer and mathematical physicist rather 
than to the pure mathematician. 
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There are now available a number of methods of evaluating opera- 
tional forms, among which may be mentioned 

The original Heaviside expansion theorem, 
Operational Division which gives a series solution, 
Contour Integration of the Bromwich-Fourier Integral, 
Carson's Integral Equation. 

It is thought that this extension to the expansion theorem will be of 
value as another way of evaluating in closed form certain operational 
expressions, especially those involving fractional exponents. 

In preparing this paper, the author wishes to acknowledge his 
indebtedness to Mr. R. M. Foster for the contribution of Equivalent 
No. 7 in its present form and for notes regarding the evaluation of the 
error function. He is also indebted to Mr. J. R. Carson for reading a 
draft of the manuscript and for a number of helpful suggestions. 



